Effects of localization and amplification on distribution of intensity transmitted 

through random media 
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We numerically study the statistical distribution of intensity transmitted through quasi-one di- 
mensional random media by varying the dimensionless conductance g and the amount of absorption 
or gain. Markedly non-Rayleigh distribution is found to be well fitted by the analytical formula of 
Nieuwenhuizen et al, Phys. Rev. Lett. 74, 2674 (1995) with a single parameter g . We show that in 
the passive random system g' is uniquely related to g, while in amplifying/ absorbing random media 
g' also depends on gain/absorption coefficient. 



Light transport in a quasi-one dimensional system is 
described by three directly measurable quantities: T a b - 
transmission coefficient from an incoming channel a to 
an outgoing channel 6; T a = ^ T ab - total transmission 

b 

coefficient from channel a to all outgoing channels; g — 
^2 T a i, - transmittanceQ,0- A wave transmitted through 

a, b 

a random medium is a coherent sum of a large number 
of contributions due to different propagation paths. In 
the diffusive regime, (g) = g ^> 1, the contributions T a b 
are largely uncorrelated. Gaussian distribution of the 
field components leads to renowned Rayleigh distribution 
P(sab) = exp (—s a b) of the normalized intensity 3], s ao = 
Tab I (Tab), with a simple relation between the moments 
of the distribution: (s™ b ) = n\(s a b) n . Such factorization 
cannot be exact 0, since it implies the complete absence 
of nonlocal correlations [7[. Kogan et al [8| related the 
distribution of the normalized total transmission s a — 
T a /(T a ) to that of s ab : 
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P(s ab ) = / —P(s a ) exp 

J s a 
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(1) 



The Gaussian distribution of P(s a ) for g 3> 1 leads to 
the deviation from Rayleigh distribution at large inten- 
sity s a b 3> 5- The deviation has the same origin as the 
universal conductance fluctuations and nonlocal inten- 
sity correlations in transport Q, Q- The distribution of 
the intensity was obtained in a closed form in Refs. I^Hfl 



P(Sa) 
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exp [xs a - $ (g',x)], 



(2) 



where $o(g',x) = g' In 2 (yi + x/g' + VW) . Eq. © 

was derived under the assumption of g ^> 1, thus g' = g. 
The expressions for P(s a ) and P(s a b) have been verified 
in experiments 0, 0, Il2l Il3j . Unexpectedly, the exper- 
iments demonstrated that Eq. (J2J worked well even for 
moderate values of g ~ 10, and in the presence of signifi- 
cant absorption. Moreover, based on the statistics of the 



transmitted intensity, the localization criterion, 

g 1 = 2/3var(s a ) = 4/3[var(s afc ) - 1] (3) 

equals to unity, was surmised|llj. Such definition of g' 
can be used, irrespective if the Eq.Q holds. However, 
if Eq.© is applicable, g' obtained from Eq.© should 
match the one obtained from the fit of the entire dis- 
tribution of s a b with Eas. (|ll2|l Jjj. O ne question to be 
addressed in this paper is whether the Eqs.(l,2) still hold 
in the regime of incipient photon localization g ~ 1, and 
if so, how the intensity distribution fitting parameter g' 
is related to the system properties. 

We will also investigate the effects of optical amplifi- 
cation and absorption on the transmitted intensity dis- 
tribution. The effect of the amplification on purely ID 
transport was studied by a number of groups 14] . In a 
quasi- ID system, the probability of reflectance by a ran- 
dom amplifying medium was also obtained by Beenakker 
et al\l^ . Zyuzin showed that the fluctuation of T a grows 
faster than its average value in a amplifying random 
medium near the lasing threshold However, it is not 
clear whether the statistical distribution of transmission 
coefficient of an active system differs qualitatively from 
that of a passive system. Brouwer demonstrated that 
i is inapplicable in the limit of strong absorption 
whereas we are interested in the weakly absorb- 
ing systems with V Dt > L, where D is the diffusion 
coefficient, and r is the absorption time. 

In this paper, we numerically calculate the statistical 
distribution of transmitted intensity P(s ao ) in a quasi-lD 
system. In particular, we check whether P(s a b) still satis- 
fies Eqs.(l,2). In real experiment or numerical study, the 
local intensity distribution P(s ab ) is obtained by collect- 
ing the data of transmitted intensity from many random 
configurations. Among them, there exist rare configura- 
tions that could lase in the presence of gain. Light inten- 
sity would diverge if gain saturation is neglected. In the 
diffusive regime g> 1, this problem is limited only to the 
immediate vicinity of the diffusive lasing threshold^]. 
For the systems we considered, g < 10. Strong fluctu- 
ation of lasing threshold results in a non-negligible per- 
centage of realizations whose lasing threshold is met even 
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FIG. 1: Numerically calculated P(s a b) for four samples. 
Dashed lines represent the fit with Eas. l|ll2^ . The inset shows 
P(sab) calculated from Bgs. lll2H with the kernel &o(g',x) 
(solid line), or <&i(g',x) (dashed line), g = 3. 

at moderate gain (i.e., not very close to the diffusive las- 
ing threshold). Although gain saturation could prevent 
the divergence of laser intensity, the actual value of the 
saturated intensity depends on the properties of the gain 
material. In order to eliminate any material-dependent 
effect on P(s a b), we disregard the contributions of the las- 
ing configurations to the intensity distribution^]. Our 
numerical studies show that the resultant P(s a b) can be 
well described by Eqs.(l,2) with a reduced g' . 

We numerically obtain P(s a b) as the histogram of the 
intensity transmitted through a quasi- ID system - a 2D 
metallic waveguide filled with circular dielectric scatter- 
ers. We employ the method developed earlier in Ref. 
Eil for passive and Ref. for amplifying and absorb- 
ing systems. In quasi-lD geometry the transition from 
diffusion g ^> 1 to localization g < 1 can be achieved 
even in the weak scattering regime by increasing the sys- 
tem length L beyond the localization length £. First, we 
study the intensity distribution in passive systems as g 
approaches 1. The random system is characterized by 
the Thouless number 6 = Sv/Av. The average mode 
linewidth 8v is obtained from the width of |Ce(Az/)| 2 di- 
vided by a numerical factor 1.46 (Ce{Av) is the spectral 
field correlation function). The average mode spacing is 

Ais = cj [ttL 1 ' Nn 2 e f j^j , where V = L + 2zb is the effective 
length of the random media, Zb accounts for the bound- 
ary effects. N is the number of waveguide modes, c is 
the speed of light in vacuum, and n e f / is the effective re- 
fractive index. All these parameters can be determined 
independently0,|2(]]. In passive diffusive system 5 coin- 
cides with the dimensionless conductance g. We consider 
four samples labeled from 1 to 4. The lengths of the 
first two samples are related as L2 — 2Li, and scatterer 
density is chosen in such a way that the increase of L 
would be offset by the change of the transport parame- 
ters, mostly by the transport mean free path I, to give 
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FIG. 2: Symbols depict 1/g' versus l/g for samples #1,3 — 
5. In the zeroth order of l/g , g'{go) = g (dotted line). 
The dash-dotted line is g'(g ) with the first-order correction. 
Dashed line plots the exact solution g(g ) from Ref. I22I The 
inset shows the linear dependence of 5 on g . The solid line 
has a slope of one. 

the identical value of 5 = 4.4. The samples #3 and #4 
had the same scatterer densities as in #1, but L is in- 
creased yielding S = 2.2 and 5 = 1.13 respectively. Fig. 
H shows that samples #1 and #2 are fitted by Eas. (|ll2|) 
with the same g' — 2.9. This means that the distribution 
of local intensity is determined by single parameter g' , 
which depends only on S (or g, see below). The samples 
#3 and #4 are fitted well with g' = 1.25 and g' = 0.4. 
Such agreement confirms the applicability of Eas. pi2p 
down to g' = 0.4, even smaller than that reported in mi- 
crowave experiments of Refs. 0, El El In what follows 
we give an argument that may shed some light on this 
agreement. 

In Ref. |^] van Langen et al found that the proper per- 
turbation series should be constructed for &(g',x), the 
kernel in Eq.JSJ). The first order (in 1 /a) correction to 
$o0',af) of Refs. H, E3 was calculated 211. In the in- 
set of Fig. n we plotted P(s a b) calculated within both 
approximations. We see that Eq. J2J) with $0 starts to 
fail for s a b < 1 region, while the asymptotic behavior for 
s ab ^ 1 is preserved. In the transition from diffusion to 
localization threshold (g ~ 1), P(s a b > 1) should be well 
described by Eas. H1I2|I with decreasing g 1 - quantitative 
change. But for g<l, P{s a b) changes to the lognormal 
distribution - qualitative change, not captured by $0 ap- 
proximation. The numerically obtained intensity distri- 
bution, indeed, shows some deviation from Eas. (|ll2| at 
small values of s a b- The deviation, however, is noticeable 
only when P(s a b) is plotted on the linear scale. On the 
logarithmic scale, the distribution has a pronounced tail 
at large s a b, where the correct asymptotic prevails. 

Next, we would like to analyze the value of g'. The 
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FIG. 3: Numerically calculated P(s a ),) for sample #3 with 
and without gain /absorption. The arrows shows the direction 
of increasing gain/absorption. Dashed lines represent the fit 
with Eqs. 11121 with g' equal to 1.5 (absorption), 1.25 (pas- 
sive), and 0.65 (gain). 



following four quantities should coincide in a diffusive 
system: g', g, 8, and go = (n/2)n e f fNl/ L' . Fig. |2Jshows 
the relations of these four quantities in the regime of in- 
cipient localization. The numerical data included sample 
#5 of 5 — 7.6. The inset of Fig. |2]shows that for the stud- 
ied systems 8 ~ g a , and thus our data for S(g ) should 
lie along the dotted line of slope 1 in the main plot. The 
rigorous solution g{g ) obtained by Mirlin in Ref. is 
shown in Fig. [5] with the dashed line. Our data for 1/g' 
lies above the dotted and dashed lines, revealing the ten- 
dency of g' to decrease below g a , <5, and g in the regime 
g ~ 1. The first order (in l/g ) correction to g'(g ) — g D 
can be obtained from the variance of the intensity distri- 
bution, var(s fl f ) ) ~ 1 + 4/3.g + 8/15g% upto the second 
order of 1 j g Q 2] . The approximation of g'{g ) obtained by 
substitution of the above expression into Eq. © is plot- 
ted in Fig. [21 with the dash-dotted line. The fact that 
the data points of 1/g' is above the dash-dotted line indi- 
cates the higher-order terms are not negligible and their 
contributions lead to a further reduction of g' . 

Next, we study the effect of coherent amplification on 
intensity distribution. Fig. [3] compares the local inten- 
sity distributions found in sample #3 with and without 
gain/absorption. The value of the absorption time is the 
same as that of the amplification time. The latter was 
equal to t = 5t c '', where r cr is the critical amplification 
time corresponding to the diffusive lasing threshold|23T|. 
Even at such low level of gain, some of the random real- 
izations lased, due to strong fluctuation of lasing thresh- 
old. The numerical results presented in Fig. [3] contain 
only the contributions from non-lasing realizations. 

First of all, we would like to point out that Eas.HH2|) 
give a good fit to the numeri cally obtained P(s a b) with 
gain/absorption. In Refs. 0, 13 it was also found that 
Eas. l|ll2|) describe well the intensity distribution even for 
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FIG. 4: The numerical data of 1/g' versus 1/8 for samples 
#1 (squares) , #3 (triangles) , #5 (circles). Solid (open) sym- 
bols represent the systems with gain (absorption). The pen- 
tagons correspond to the passive system. The dotted, dashed, 
and dash-dotted lines are the same as in Fig|2] Solid lines 
represents g'(S). The inset shows the explicit dependence of 
g' on the amplification/absorption time r normalized to r cr . 
Dashed (solid) line plots Eq.JlJ with positive (negative) r. 



absorbing systems. Here, we present a systematic study 
of the effect of absorption on P(s a b) and also consider, 
for the first time, the effect of amplification. Fig. |3| 
shows that the presence of gain leads to an increase of 
P(sab) in the regions s a b *C I and s a b ^ 1 and, therefore, 
an enhancement of intensity fluctuations. The effect of 
absorption is exactly the opposite. From this qualitative 
analysis and Eq. , we see that amplification reduces the 
value of g' , while absorption increases it. 

Fig. 01 shows the effect of amplification/ absorption 
on g' for samples #1,3,5. We can still define 5 for 
amplifying/absorbing systems through the correlation 
functions [19(. In contrast to g a , 8 depends on the amount 
of gain/absorption. The change of 8 is directly related 
to the change of the correlation linewidth 8v, while Az/ 
remains the same as in passive system. The superlinear 
increase of 1/g' with 1/8 in Fig. ^demonstrates that g' 
decreases significantly faster than 8 with increasing gain. 
This indicates that intensity fluctuations are more sen- 
sitive to amplification than the average mode linewidth 
8v. The long propagation paths extract more gain, lead- 
ing to an enhancement of non-local correlations 19] . The 
enhanced correlation results in a further deviation from 
Rayleigh statistics. The absorption, on the other hand, 
causes a reduction of fluctuations. The sublinear decrease 
of 1/g' with 1/8 in Fig. 0] reveals that g' increases slower 
than 8 with increasing absorption. The dependence of 
g' on the absorption time r can be determined pertuba- 
tively from the known result for var(s a h) in an absorbing 
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TABLE I: Comparison of localization threshold criteria in 
passive system. 

jj_ 9o 5 g_ 

1.0 2.0 2.0 1.3 

0.3 1.0 1.0 0.43 

0.7 1.6 1.6 1.0 



but increases sublinearly with absorption constant. 
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system of g a » 1 
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3ff'(ffo,i") 



var(s afc ) - 1 = - — A 2 {t) - 

3g 



Mr), (4) 



where A 2 {t) and A 3 (t) are given in Refs. [13, |2J. The 
expressions for Aj were derived for absorbing systems, 
where r is positive. In the negative absorption model, r 
becomes negative for systems with gain. In Ref. the 
applicability of the negative-r expressions was discussed. 
Obviously, the contribution of the lasing realizations 
should be omitted to avoid the divergence of var(s a h). 
For the samples #1 and 5 that have larger g a , we com- 
pared our data to Eq.QJ in the inset of Fig. 0] Also, 
5{g ,T) can be found from the width of |Cg(Az/)| 2 ^|. 
By eliminating r from g'{g ,T) and 5{g ,r) we obtained 
g'(5), shown as the solid lines in Fig. 0] The deviations 
of our data points from the solid line increase with the 
decrease of g , because the contributions of the higher 
order terms in l/g cannot be neglected. 

In Ref. ^3 g' was proposed as the localization crite- 
rion parameter. Table ^ shows that in passive systems 
the criteria based on g' = 1, g a = 1, g = 1, or S = 1 
differ by only a numerical factor of order one. In the 
presence of gain/absorption the situation is essentially 
different. By definition, g a only contains information of 
passive system, whereas g' and S account for the effect 
of amplification/absorption. In this case the difference 
between g', S, and g, is not merely numerical. Am- 
plification/absorption results in a decrease/increase of 
g' and S, but an increase/decrease of g. With the in- 
crease of absorption, S increases without a bound, while 
g'(r — > 0) — * Ag'{r — > oo)/3 in the limit g > 1. In 
sharp contrast, we see that in amplifying systems g' di- 
minishes superlinearly with r cr /r (inset of Fig. 0J, while 
S decreases almost linearly with r cr /r as reported in Ref. 
Hl^ Our numerical result in the inset of Fig. Q]also sug- 
gests that g' fall below unity prior to <5 in an amplifying 
system. Therefore, g' is more sensitive to amplification 
but less sensitive to absorption than S. 

In conclusion, we numerically calculated the statistical 
distribution of transmitted intensity in quasi-lD random 
medium. In a passive system P(s a b) is well described 
by Eqs.jnU down to g' = 0.4, far beyond the g' > 1 
regime where it was originally derived. Surprisingly, Eqs. 
(HIL'I) also hold for amplifying/absorbing systems with g' 
different from the value of the passive system. Our data 
show that g' decreases superlinearly with gain constant, 
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